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Abstract 

The influences of quantum effects on the structure of the phase-space of 
generalized Gauss-Bonnet theory, introduced by the Lagrangian F(R,G), 
have been studied. G is the Gauss-Bonnet invariant, and the quantum 
effects are described via the account of conformal anomaly. It has been 
shown that the quantum effects change many aspects of the attractors 
of F(R, G) gravity models in the R ~ H plane, including the location of 
the attractors, the number of them and their stability properties. These 
variations are not, in general, from the type of small perturbations, but 
instead, it can induce the great, so not ignorable, variations which have 
root in the "singular perturbation" nature of this effect. In other words, 
one can not ignore the quantum corrections and must be always consid- 
ered. The influences of the perfect barotropic fluids on this problem have 
been studied, and it has been shown that this kind of matters do not alter 
the quantum effects. It has been shown that the classical contribution 
of the coupled-quintessence model, which is responsible for inducing the 
quantum effects, is of this type, that is a perfect barotropic fluid, and 
therefore can not change our results. 

1 Introduction 

Independent observational data indicate that we are now in accelerating phase 
of the universe. The supernova observations directly support this accelerating 
expansion, and the microwave background, large scale structure and its dy- 
namics, weak lensing and baryon oscillation observations indirectly verify this 
phase [1]. It is believed that nearly 70% of our present universe is composed of 
dark energy, the physical object which is responsible for the effective negative 
pressure required for accelerating cosmic expansion. 

There are two main dynamical models that explain some features of dark 
energy. The first one is based on the standard Einstein theory, when some other 
physical terms added to it. The simple cosmological constant model, with no 
dynamics [2], the scalar field models, including quintessence, phantom, quin- 
tom and hessence models [3], and some other models like the k-essence models, 
chaplygin gas models, and holographic dark energy models [4] are examples of 
the first class of dark energy models. 
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The second class of dark energy models are those based on the assumption 
that the gravity is being (nowadays) modified. These modified gravity theories 
can be divided to two main subclasses. The first subclass consists of the models 
in which, besides the scalar-curvature R, there exists a scalar field (or either a 
vector field) which interacts nonminimally to gravity. These models arc known 
as the scalar-tensor models, where after the simplest Brans-Dicke theory [5], 
several generalizations have been introduced [6]. Recently it has been shown 
that the generalized scalar-tensor theories admit the co = — 1 and deceleration 
to acceleration transitions, the phenomenon which can be affected by quantum 
effects [7]. The equation of state parameter to is defined through to =p/p, where 
p is the pressure and p is the energy density of dark energy. In cosmological 
constant model, to takes the constant value lua = — 1, and in dynamical models, 
it is a function of time, i.e. oj = to{t). 

The second subclass of modified gravity models are those whose actions are 
not the simple term R, but instead they are, in general, an arbitrary function of 
all algebraic invariants built up with the Ricmann tensor, including the scalar 
curvature R, the quadratic invariants P = R^R^ and Q = R 01 ^ 6 R a ^g, 
and other independent invariants of higher orders. The f(R) gravity, whose 
Lagrangian is an arbitrary function f(R), is the simplest and the most famous 
modified gravities of this subclass, i.e. 



S = (fx 



(1) 



in^ = c = G= l units, k 2 = 8ir and C m is the Lagrangian density of dust-like 
matter. Many aspects of f(R) gravity, including their local properties, have 
been studied. See [8] and references therein. 

Restricting ourselves to quadratic invariants, the generalized modified grav- 
ity models are defined through 

S = J d A xV=g [f(R,P,Q)+C m ], (2) 

where f(R,P,Q) is an arbitrary function. Some aspects of these models, in- 
cluding their attractor solutions and their stabilities have been discussed in [9] . 
Studying the propagators of these models has shown an important result, i.e. 
the generalized modified gravity models have, in general, the graviton ghost, 
unless it satisfies [10] 

(/p + 4/ q )U=h =0, (3) 
where the subscripts denote the partial derivatives, e.g. 

_ df(R,P,Q) 

IP - gp , (4) 

and the curvature Rq is defined by the equation 

/ - \RqSr - \Ro 2 fp - -V/q = 0. (5) 

The above relation is the equation of motion of the metric at constant scalar 
curvature R = Rq. In eq.®, Q — i?o 2 /4 and P = Rq 2 /6. Now if one considers 
the Gauss-Bonnet (GB) invariant G 

G = R 2 -AP + Q = R 2 - AR^R^ + R^R^. (6) 
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it is easily seen that the condition ([3]) satisfies, and this is a reason why the 
generalized GB gravity, defined by 



is an important candidate of modified gravity theories. 

The generalized GB gravity, where the f(R) gravity ([T]) is a special example 
of it, has been first introduced in [11]. The hierarchy problem of particle physics 
and the late time cosmology in the context of F(R, G), have been studied in [11] 
and [12], and its behavior under phantom-divide-line crossing and deceleration 
to acceleration transitions has been investigated in [13]. Recently, the two- 
dimensional phase-space of the generalized GB models, i.e. the R — H space 
where H is the Hubble parameter, has been studied in [14] and the various 
aspects of the de-Sitter attractors of these models have been discussed. 

The present paper is devoted to study the contribution of quantum effects 
to the attractors of generalized GB gravity. The quantum effects arc described 
via the account of conformal anomaly, reminding about anomaly-driven infla- 
tion [15]. The contribution of conformal anomaly in energy conditions and 
big rip of phantom models has been discussed in [16], and its influence on the 
u> = — 1 crossing and deceleration to acceleration transition of quintessence and 
phantom models, the F(R, G) gravity and the generalized scalar tensor models 
have been discussed in [17], [13] and [7], respectively. 

The late-time behavior of dynamical models, including the dark energy mod- 
els, is an important problem, both in mathematics and physics, and is studied 
in a framework known as the attractor solutions of dynamical systems. Many 
properties of attractor solutions of dark energy models have been studied [18]. 
Two important parts of attractor studies of any dynamical system are: a) choos- 
ing the phase-space of the model and b) obtaining the location of attractors and 
specifying the stability of them. In F(R, G) gravity models, the only possible 
choice of phase-space for general F(R, G) model is R — H space. So attractors, 
with properties R — and H — 0, or R = R c and H = H c {R c and H c are 
constant values ), are de-Sitter solutions of generalized GB models. Now, as we 
will show, the quantum effects have important, and non-ignorable, contributions 
to this problem. It changes the position of the attractors and their stability be- 
haviors, and, for some cases, it produces the new attractors where some of them 
do not lead to classical attractors, a phenomena known as "singular perturba- 
tion" in mathematics. This shows that the quantum effects are not ignorable, 
i.e. they do not only lead to small corrections for classical solutions. For some 
other cases, it can remove the degeneracy of solutions, which is known as "bi- 
furcation" in the mathematics of dynamical systems. And, interestingly, it can 
transform the critical curves, an infinite number of stable attractors locating on 
R = 12H 2 curve in the phase-space of some specific case of F(R, G) models, to 
a unique point. 

The scheme of the paper is as follows: In section 2, we discuss the crit- 
ical points and the stability conditions of F(R, G) models in the presence of 
quantum correction terms. It is shown that in the H — ► limit, the relations 
are correctly reduced to classical relations. In section 3, we discuss the various 
examples of f(R) and F(R, G) models. The new aspects of the phase-space of 
these models, including those introduced in the last paragraph, can be seen via 
these examples. It is shown that the numerical calculations verify our results. 
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The usual matter's contributions to the location and behavior of attractors are 
discussed in section 4. It is shown that the barotropic fluids can not change 
the structure of the phase-space of F(R, G) models. It is also shown that the 
classical contributions of the coupled-quintessence model, which is responsible 
for producing the quantum correction terms considered in section 2, can not 
alter our results, which is an important consequence. 



2 Quantum corrected critical points of F(R, G) 
gravity 

Consider the generalized GB dark energy with action Q. Variation of this 
action with respect to the metric g^ results in [11] 

hp^ + ^ g ^F(R, G) ~ 2F G (R, G)RR pv + 4F G {R, G)R" p R lJp 

- 2F G (R, G)R waT R v paT - 4F G {R, G)R PP ™ R pa + 2(W' 1 WF G (R, G))R 

- 2g p »(V 2 F G (R, G))R - 4(V p V"F G (i?, G))R" P - 4(V p V l/ J F G (i?, G))R PP 

+ 4(V 2 F G (i?, G))R p,J + 4g pu (V p V a F G (R, G))R pa - 4(V p V CT F G (i?, G))R ppva 

- F R (R, G)R pu + V^V V F R {R : G) - g^V 2 F R {R, G) = 0. 

(8) 



T pv is the energy-momentum tensor, defined through 

2_6S 1 



Tfiv = (9) 



and the subscripts of F(R, G) denote the partial derivatives. For the back- 
ground metric, we consider the flat Friedmann-Robertson- Walker (FRW) met- 
ric, defined by 

ds 2 = ~dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). (10) 

(t, x, y, z) are comoving coordinates and a(t) is the scale factor. The three 
independent Friedmann equations of F(R, G) gravity, then become [13] 

-6H 2 F R (R, G) = F(R, G) - RF R {R, G) + QHF R (R, G) 

+ 2AH 3 F G (R,G)~GF G (R,G)- Pm , (11) 

R = 6(H + 2H 2 ), (12) 
G = 2AH 2 {H + H 2 ), (13) 
and the evolution equation of matter field is 

p m + 3H(p m +p m ) =0. (14) 

H = d(t)/a{t) and "dot" denotes the time derivative. 

To study the quantum effects, we consider the following standard Lagrangian 
of a scalar field 4>, known as the coupled-quintessence model, in which <j) is 
nonminimally coupled to gravity 

C = \ (- - W) 2 ) , (15) 
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where the renormalizibility requirements force £ to be £ = 1/6. Calculating the 
effective action of this conformal invariant Lagrangian at one-loop level, results 
in a nonvanishing trace for the energy-momentum tensor, which is classically 
traceless. This trace, i.e. the trace/conformal anomaly, is [15,19] 

T A = b(F + + b'G + b"UR . (16) 

The subscript "A" denotes "anomalous", G and R are Gauss-Bonnet and Ricci 
scalars, and F is the square of the 4d Weyl tensor 

F=\R 2 - 2R flv R» v + R^ al3 R^ al3 . (17) 

b, b' and b" are given by 

_ 3(iV + 6N 1/2 + 12JV a + eiliVa - 8iV HD ) 

b= ^60(4^ ' 

N+ UN 1/2 + 627V! + 14117V 2 - 28N RD { ' 

360(4tt) 2 , & - U. 

Ea. (fT8|) is for the cases where there exist N scalars, Ni/ 2 spinors, Ni vector 
fields and N 2 (= 0, 1) gravitons, and A^hd higher derivative conformal scalars. 
The energy density pa and pressure pa, corresponding to Ta, can be found by: 



1 



PA = -— / a 4 HT A dt, (19) 
a Jo 

and 

PA — ^ (T A + Pa) , (20) 
respectively. The resulting relations in FRW metric are [16] 

p A = -6b' H 4 - (h + b")(-6HH - 1%H 2 H + 3H 2 ) , 

2 (21) 
p A = b'{6H 4 + 8H 2 H) + (-b + b"){-2H - 12HH - 18H 2 H - 9H 2 ). 

Note that in SI units, the above relations have h/c 3 as the prefactor. It can be 
easily shown that the above expressions satisfy 

Pa + 3H( Pa + pa) = 0. (22) 

The natural method for computing the quantum corrections in F(R, G) gravity 
models is to add pa to Friedmann equations, i.e. to change p rn in eq. fTTT) to 

Pm + PA- 

Using eas. (Tl2|) and (fl~3|) . one finds 

G = AH 2 (R - 6iJ 2 ), (23) 

from which 

G = -HR 2 + 192H 5 - 32RH 3 + 4H 2 R. (24) 
3 
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Also ea.([2l"]) leads to 

PA = -6b' H 4 + (lb + 6") (hR + RH 2 -^) ( it, i 

The time derivatives in ea. (|ll[) can be also written in terms of R and G, using 

j t f(R,G) = f R R + f G G. (26) 
In this way, R and H are found from Friedmann equations as following 

R = ; 5 ; : To — X { (R - 6H 2 )F R + GFq 

6H{F RR + m 2 F RG + \6H±F GG )~{lb + b")H \ V ' 



- 288H 2 (- - 2ff 2 ) (F flG + AH F GG ) + (-b + b"){RH 2 - — ) - 6fo'tf 4 + p m r . 

(27) 

H=--2H 2 . (28) 
6 

These are the set of autonomous equations, where its phase-space is two dimen- 
sional (R — H) space. Here, for simplicity, we do not consider the matter field, 
i.e. p m — 0. We will come back to matter field in section 4. Note that in the 
right-hand-side of ea. (j2"7"]) . the relation (|23j) must be used for G. 
Our autonomous equations are in the form 

H = fi(R, H), 

. (29) 
R = h(R,H). 

The critical points are found by setting eas. (l2"T]) and l|28p equal to zero. The 
results are 

~RF R + GF G -F = 6b' H 4 , (30) 

and 

R=\2H 2 , (31) 

from which ea. (|23|) results in 

G = 24H 4 = — (32) 
6 v ; 

at critical points. The above equations must be solved to obtain the coordinates 
of critical points in R — H space, i.e. (R c , H c ). Note that in b' — > limit, 
which is, in fact, the h — > limit (mentioned after ea . (|2~1"T) ) . the critical point 
relations (f3T)]) and (|31[) are reduced to classical relations obtained in [14]. It is 
also interesting to note that from three quantum parameters b, b' and b" , only 
the parameter b' appears in critical point equations (|30p and (|3ip. 

To study the stability of each critical point, the eigenvalues of the matrix 

M = ( d J}l™ ) (33) 



dh/dH dh/dR 



R = Rc ,H=H c 
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must be calculated. The index "c", denotes "critical value". The critical point 
is a stable attractor if the real parts of all eigenvalues are negative. The stable 
attractors divide to two categories: it is "node" when the eigenvalues are real, 
and "spiral" if they are complex conjugate. For real eigenvalues, but with 
opposite sign, the critical point is a saddle point. 

For autonomous eqs. (|2"7]) and l|28p. the matrix M becomes 



M 



-AH 1/6 
S H 



R=R C ,H=H C 



where 



S 



-6F R + (U + b" -b')R 



3{F RR + 8H 2 F RG + 16H 4 F GG ) 



lilb+wy 



The eigenvalues of M then become 



Al 2 — — 

1-2 2 



-3H± \ (3H) 



24H 2 ) 



(34) 



(35) 



(36) 



J R=R c ,H=H a 



To ensure a stable attractor, the real parts of Ai and A2 must be negative, which 
satisfies if, and only if, S + 24H 2 < 0, or using eq.(3T]): 



S = S + 2R 



< 0. 



R=R C ,H=H 

At b' —>■ limit, the above condition, multiplied by —1/6, is reduced to 

Fr 



V 



3{F BR + 8H 2 F RG + 16H*F GG \ 



- 4H Z 



> 



(37) 



(38) 



R=R C ,H=H C 



which is the same as the classical stability condition, obtained in [14]. It must be 
noticed that in contrast to the critical point equation (|30|) which only depends 
on 6', the stability condition ([3T)l depends on all three parameters b, b' and b" . 

Before considering the complicated examples in the next section, it is seen 
that for i?-independent Lagrangian 



F(R,G)=F(G) 



(39) 



where classically has no stable attractor, since r\ = —AH 2 < 0, the quantum 
effects can induce the stable attractors for these models. For example when 
F G g — 0, ea. ([37f becomes 



2b' 



b" 



-R 



(40) 



which is negative for the most background fields, i.e. eq. (fT5|) shows that b' < 0, 
except for 287V H d > A r + lliV 1/ 2 + 62iVl + 1411A? 2 - So many of the F{G) gravity 
models may have the stable quantum attractors. 



3 Some examples of f(R) and F(R, G) gravities 

In this section we study some specific examples to explore some important fea- 
tures of quantum attractors. 
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3.1 F(R, G) = ai + a 2 G In G models 

As the first example, we consider a specific F(G) gravity model with action 
ai + a 2 G In G. The classical critical points are found by ea. (|5D|) . when b' sets to 
zero, 

GF G -F = ^ a class - = — , (41) 

(42) 

G F a F = -G : ■:■ <7,.' !m ""' ' : — — t ■ (43) 

(44) 



or, using eq. 

class. 

At quantum level, eq. (|30p results in 




from which 



quant. 



24ai 



4ai 



4a 2 



4a 2 - &' 

Depending on the parameters a±, a 2 and b', different situations may arise. For 
example if a\a 2 > 0, so that i? c class - in ea.(|42p becomes a real positive number, 
the classical critical point is an unstable attractor (see after ea.(|39p). but the 
quantum critical point can be a stable attractor. In Figs. 1-3, the results of the 
numerical calculations of eas.(|27p and (|28p for F(G) = 1 + 2GlnG model has 
been reported. Because of the numerical factors of b and b' in ea.(|18p. it is easier 
to rescale F(R, G)s by c, where c" 1 = 360(4tt) 2 . Therefore by F = 1 + 2GlnG, 
we mean F = ,a (1 + 2GlnG). Now we consider two different background 

360(47r) v j ° 

matters. We take N 1/2 = Ni = N 2 = N KD = 0. For N = 1, one has b = 3 
and b' — —1 and for N = 5, b = 15 and 6' = —5. In Fig.l, It can be seen that 
R c class - — y / 6ai/a 2 = Vo is not a stable attractor, Fig. 2 shows that in N = 1 
case, the quantum corrected attractor ^quant- 
not yet stable, but in the case N 
stable attractor. 



R 

40 : 

30 
20 
10 



^24ai/(4a 2 - V) = yf2Aj§ is 
5, Fig.3 shows that R^ uant - = ^24/13 is a 



2 4 6 8 10 
Figure 1: The plot of R(t) of F(R, G) = 1 + 2G In G model at classical level. 

Other interesting cases also exist. For aia 2 < 0, where classically there is no 
critical point, quantum effects can induce a stable attractor. Fig. 4 shows the 
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Figure 2: The plot of R(t) of F(R, G) = 1 + 2GlnG model, with N = 1. 




0.334 0/335 



Figure 3: The spiral paths in i? — H plane of F(R, G) 
N = 5, at quantum level. 



1 + 2GlnG, with 



phase-space of F(R, G) = 1 — |GlnG model with N = 1 background. As it is 
seen, the point (R C ,H C ) = (6.93,0.76) is a stable attractor of this model, which 
has no analogous at classical level. 




Figure 4: The paths of F(R, G) = 1 — |lnG, with AT = 1, at quantum level. 



3.2 F(R, G) = oi + a 2 # + 03-R 2 models 

In the previous example, there is one critical point in each classical and quantum 
mechanical levels, and at b' — > 0, i.e. H — ► 0, limit, i? c quant - — > i? c class '. Now 
we consider the class of models in which the number of classical and quantum 
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mechanical critical points are not the same, and more, at b' — > limit, some of 
the quantum attractors do not lead to classical ones. 
Consider the following f(R) gravity model 

f{R) = 01 + a 2 R + a 3 R 2 . (45) 

The only classical critical point is found by ea. (j30|) . with b' = 0, as following 

R c cUss =-—. (46) 
a 2 

The condition of stability of this solution, using eq. ([55]) . is 

V=-^>0. (47) 
6a 3 

So to have a positive curvature attractor, a\a 2 must be negative, and for its 
stability, a 2 a 3 must be positive. Note that because of eci. pTj) , the negative 
values of R c are not acceptable, since they lead to imaginary Hubble constants, 
which are not physical. 

After adding the quantum correction terms, the action (|45[) has two quantum 
critical points. Using (f5T7|) . one finds 



, -6a 2 -V(6« 2 ) 2 -24 fll 6' 

n cl ~ 



quan , _ ~6a 2 + V(6 fl2 ) 2 -24a^ 

U c2 - y 



(48) 



and their corresponding Ss, which determine the stability behavior of these 
points, are 

2 A /(6a 2 ) 2 - 24ai6' 

Si = -S 2 = -Y- -* . (49) 

12a 3 -(§& + &") 

Let us first study the b' — ► limit of curvatures in ea.(|48p. which their values 
depend on the sign of a 2 : 

^quant. _^ _12o 2 ^ ^quant. _^ _ 2fli = ^ class . ^ for ^ > Q ^ 

a 2 

and 

i? cl quant •— > -— = i? c class ' , i? c2 quant - — - , for a 2 <0. (51) 
a 2 b 

Now we encounter a new situation, one of the quantum curvatures goes to 
classical value, and the other one diverges. 

This behavior is the characteristic of " singular perturbation theory" . Math- 
ematicians divide the perturbation theory to two categories, "regular" and "sin- 
gular" . In regular perturbation problem, the solution of the perturbed equation 

/i(x) + e/a(aO = (52) 
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IS 

oo 

x(e) = ( 53 ) 

n=0 

where do, i.e. the zero-order solution, is the solution of the unperturbed equa- 
tion: 

/i(oo)=0. (54) 

This is the usual behavior that we almost always encounter in physics. But in 
some perturbation problems, the situation is different. In singular perturbation 
problem, these two solutions, i.e. the "zero-order solution" and the solution of 
"unperturbed problem", do not coincide. In fact, the zero-order solution may 
depend on e and may exist only for nonzero e. This situation occurs whenever 
the power of the perturbation terms are greater than the unperturbed terms [20] . 
Consider, for example, the following algebraic equation: 

ex 2 + x - 1 = 0. (55) 

Their solution are 




l±vT+4e) ^\ 1 ; ' ' (56) 



2e 

At e — > limit, x\ — ► 1, which is the solution of the unperturbed equation 
x — 1 = 0, while X2, which diverges at e — > 0, has no unperturbed analogous. 

This highly dependence of a problem to the perturbation is frequently en- 
countered in chaotic dynamical systems. The appearance of this kind of solu- 
tions shows that one can not discard the perturbative terms, i.e. it can not be 
ignored, sets to zero, in the equations and therefore in the solutions. 

For the example in hand, i.e. the Lagrangian (|45p , the critical point equation 
([50)1 becomes 

- ai -~a 2 R=^R 2 , (57) 

which is very similar to eq. (|55p . So it is natural we encounter the singular per- 
turbation problem, with its mentioned behaviors. In this problem, one classical 
solution, increases to two quantum mechanical ones, which one of them, de- 
pending on the sign of 02, goes to classical solution, and the other one blows 
up. This shows that one can not ignore the quantum correction terms and they 
must be always added to classical equations of motion. The classical attractor 
can be stable or not, depending on the sign of 0.2013 , see eci.([4"T|). but in quantum 
case, one of the solutions is always stable and the other is unstable, see eq.(gSJ). 

We may have, depend on the parameters, one positive R c class - and one posi- 
tive i? c quant - , no positive R c class - and one positive i? c quant - , or one positive R c class - 
and two positive R c quant \ As a specific example for the latter case, we consider 
the following action 

F(R) = -l + R-R 2 (58) 
in which we again write in c - 1 = 360(4tt) unit. Using eqs.([36 > (|gg ]) . one finds 

^class. = 2 with n = _ ly / 6) R quant. = g + ^ ^ = _ 1 /j2/ 7) and 

i? c2 quant - = 6 - VU with 6 2 = VT2/7. So i? c class - and i? c2 quant - are unstable 
critical points and tend together at b' — > limit, while i? cl quant ' is a stable 
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attractor which blows up at this limit, see eg. ([50)) . We have chosen N = 1, and 
other iVj's equal to zero. Fig. 5 shows the behavior of R(t) in classical regime, 
which verifies the unstability nature of i? c class - , and Fig. 6 shows that i? cl quant ' 
is a stable attractor. 




5 10 15 20 
Figure 5: The plot of R(t) of F{R) = -1 + R - R 2 model at classical level. 




Figure 6: The quantum attractor of F(R) = —1 + R — R 2 model, when N = 1. 
The stable attractor R c \ = 6 + Vl2 blows up at b' — * limit. 



3.3 F(i?, G) = ai + a 2 -R + «3-R 2 - a^R 2 In i? models 

As other interesting example, we consider the following f(R) action: 

f(R) =ai +a 2 R + a 3 R 2 - a 4 R 2 In R. (59) 

This model has two classical and two quantum-mechanical critical points, with 
coordinates 



R cl 



class. 



— a 2 — y (a 2 ) — 8aia4 



2a4 



R 



class. 



-0,2 



+ \J (a 2 ) 2 - 80104 



(60) 



<;2 



2a4 
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and 



quan , - tt 2-V( a 2) 2 -8a 1 (a 4 + 6Vl2) 
- - 20^/6 (61) 

fl qua „t. _ + y/faf-Sa^ + V/U) 



L c2 



2a 4 + &'/6 



respectively. As 6' tends to zero, irrespective of other parameters, it is clear 
from eqs. (|6"0j) and (|6ip that 

i? CJ quant ' b '^° . i? CJ class - , i= 1,2. (62) 

However, this occurs whenever both classical and quantum critical curvatures, 
separately for each i in eq. (|62[) . belong to the same region, i.e. positive or 
negative curvatures. But if they are placed in different regions, e.g. i? cl class - < 
and i? cl quant ' > 0, we encounter a new problem. As a specific example, if we 
demand 

Rc desB, < q ; R c2 class - > , R cl quant - > , i? c2 quant - > 0, (63) 
then the parameters must satisfy 

6' 3a 2 

o'<0 , a 2 >0 , 0<a 4 <-- , 24^T2^ ai< °- (64) 

Because of the condition < a 4 < —b'/12, at 6' — > limit, the parameter a 4 
must also goes to zero. So the denominators of i? cl class - and i? cl quant ' both go 
to zero, with different signs, while their numerators remain fmitcQ Therefore 



class. 



i? cl quant ' b '"° » +oo (65) 

This is another interesting behaviors of quantum F(R, G) gravity models. The 
quantum effects increase the number of positive attractors, but the excess at- 
tractor blows up at h — > limit and does not approach to their corresponding 
classical solution. As an explicit example belongs to this category, we consider 
the following Lagrangian 

f(R) = -1 + 2R+ ^R 2 - ^R 2 \nR , (66) 

and choose N — 3 and other iVj's are zero, so that b' = —3. The critical points 
then become 



(Rcl 



i? cl class ,i? c2 class ) = (-3(2 + 4/V3), 3(-2 + 4/V3)j , 

! " ant -,i? c2 quant ) = (2(6 + V30), 2(6-V30)V (67) 



So we have one acceptable (positive) classical and two quantum mechanical 
critical points. All the critical points are stable attractors. i? c2 quant ' tends to 
R c2 class - at b' — > 0, while i? cl quant ' has no classical analogous. Figs. 7 and 8 show 
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Figure 8: The quantum attractors of f(R) = -1 + 2R+ |i? 2 - \R 2 In R model, 
when N = 3. 



the behaviors of the stable attractors of this model at classical and quantum 
regimes, respectively. 

The other interesting case of the Lagrangian ([59")) , is one in which the classical 
curvatures become degenerate. This occurs when (02) 2 = Ba^a^. In this case, 
one has 

D class. p class. r> Q 2 /fiQ\ 

H cl — U c2 - -"class. - ) i bS i 

and 

t-> quant. -^class. quant. -^class. 

R cl H = _ , i? c2 q = _ . (69) 

1 _ Y _ 12oI ^ + Y — 12al 

This phenomenon, in mathematics, is known as "bifurcation". In dynamical 
systems, if the problem depends on one (or more) parameter(s), and by contin- 
uous varying the parameter(s), two of the critical points collide each other, it is 
said there is a local bifurcation. In our problem, b' plays the role of the dynam- 



^^Note that for _R c2 cIass - and i? c2 quant ', also the denominators go to zero, but at the same 
time, the numerators tend to zero, such that i? c2 class ' and i? c2 quaIlt ' remain finite. 
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ical parameter. For b' ^ 0, we have two attractors By decreasing b' and 

reaching 6' = 0, they collide and a single classical attractor is produced. 

Another characteristic of this degenerate solution is its ?7-value. To study 
the stability behaviors of i? c lass.> one must calculate r\ in ea. ([3"5)) . For f(R) in 
ea. (l59l) . r\ becomes, using ea. (j3"Tj) . 



_ }r _ R _ a 2 + 2a 4 i? 

V ~ SfiiR 3 ~ 3(2a 3 - 3a 4 - 2a 4 IniZ) lR=Rc 1 j 

So for i?ciass. in cq. (|68[) . one finds t]\r=r c1b „ s =o- This means that in this case, the 
eigenvalues of matrix M are Ai = 0, X2 = —3H (see eq. (|3l)l) for classical case, in 
which S + 24H 2 must be replaced by — 677) . The appearance of zero eigenvalues 
in dynamical systems is an important point, which one of its consequences is that 
we can not predict the stability behavior of the attractors. In this case, the first 
order variation, which results in the stability matrix (|33p . is not adequate and 
one must consider the higher order approximations. This subject, for G) 
gravity models, has been discussed in [14]. 

As an explicit example of the bifurcation point, we consider the following 
example: 

f(R) = 2-4R + 2R 2 - R 2 \nR (71) 

and N = 1. Then i? class . = 2 and Rf ant - = (24 ± 2a/12)/H. i? c i ass . and Rf ant - 
are unstable attractors, while Jj^ uant - i s a stable attractor, see Fig. 9. 




Figure 9: The quantum critical point i^ uant - f f(R) = 2 - 4R+ 2R 2 -R 2 lnR 
model with N = 1. 

As pointed out earlier, because of the degenerate nature of the classical 
attractors of this example, one can not easily determine the stability behavior 
of these critical points. In fact, if one follows the method described in [14], it 
can be shown that the next to leading order variation of Lagrangian (|71|) , near 

^class. — 2, is 

U = 6.9U 2 + ■■■ , (72) 

where the positive coefficient of W 2 -term, proves the unstability of classical at- 
tractors near i? c iass. = 2. 
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4 The matter field's contributions to attractors 



In this section we want to study the contribution of matter fields to the structure 
of the phase-space, and investigate whether our results in previous sections are 
changed. 

By considering p m , the autonomous system of equations then consists of 
three equations, i.e. eas. ([27l) . (128)) and ea. (fT4|) . and therefore the phase-space 
becomes 3-dimensional, i.e. (R, H, p rn ). It must be noted that the variable p m in 
cq. CB)) can be expressed, by using the equation of state of matter p m — f(p m ), 
in terms of p m , so it is not an independent dynamical variable. 

We restrict ourselves to barotropic fluid in which 

Pm = u m p m , (73) 

and w m is a constant. For example for dust, u>d — 0, and for radiation, lo t = 1/3. 
In this case, the third autonomous differential equation becomes 

=-3H(l + u m )p m . (74) 

One can write the new system of autonomous equations as: 

H = h(R,H,p rn ) , R = f 2 {R,H,p m ) , p m = f 3 (R,H,p m ). (75) 

The critical points are found by setting H = R — p rn — 0. In this way, the 
critical point equations become the same as previous ones, i.e. eas. (|30p and 
(|3ip. and the following extra relation: 

Pmc = 0. (76) 

So the presence of the matter fields does not affect the position of the attractors, 
which comes from eas. ([30l) and (l3T|) . 

The stability matrix M now becomes a 3 x 3 matrix, as follows 




M= | S H B ) , (77) 

R=R c ,H=H c ,Pm=p m c 



in which 

B = 



1 



6H(F RR + 8H 2 F RG + 16H 4 F GG ) - (§6 + b")H 



(78) 



R — R. c .H — H c 



The eigenvalues are: Ai and A2 in eq. (|36p . and 

A 3 = -3(1 + uj m )H c . (79) 

So, as long as 

w m > -1 , (80) 

the eigenvalue A3 is negative and the stability situation is completely determined 
by eci. p7p . or (|38p . In this way we obtain an important result: For barotropic 
fluids with io m > —1, the structure of the phase-space is not changed. Note that 
the condition (|80p satisfies by all ordinary, i.e. baryonic, matters. 
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The next point that we must consider is that the Lagrangian (| 1 5[) has two 
kinds of effects in our problem. One of its role is its quantum contribution to 
the F(R, G) gravity, which we have considered it by adding its pa and pa to 
the Friedmann equations. The second role of (fT5|) is its classical contribution to 
our problem, which we have not yet considered it. To do so, one must consider 
the Friedmann equations of combined actions ((?]) and (fllj)) . i.e. 



X\ 



S — jd~ , , . , ~ , , , t 



F(R, G) + C m + i {-£Rtf - (V0) 2 



(81) 



The other way is to consider the effects of extra action (|T5|) , by its energy density 
and pressure, i.e. and p<f,, as the source terms for Friedmann equations (|2~T|) 
and (|28]) . This can be done by considering the Friedmann equations of the 
scalar-tensor model (|15[) in the background of ordinary Einstein model, i.e. 

S = \ j d^x^—g [R - iRtf - (V0) 2 ] , (82) 

which are [7] 

3H 2 = P4 , = 3£,H 2 (f) 2 + i</» 2 + 6il£# , (83) 

- 2H = P<t> + n = -2£H<j) 2 + 2iH^ - 2£# + (1 - 2£)0 2 , (84) 

and 

4> + 3H4> + £(j)(6H + 12H 2 ) = . (85) 

is therefore given by eq. (|8"3")) . p<p can be found by subtracting (|53")) from (fgl)) . 
using (|85[) . which results in 

Pl p = [2£(6£ - 1)H + 3£(8£ - l)ff 2 ] (j) 2 + 2£ff# + (i - 2^ 2 . (86) 
Then the trace of energy- momentum tensor, T^, becomes 

T = - P4 , + 3 P0 = (1 - 6£)0 2 + 6^(6^ - 1){H + 2H 2 )tf. (87) 
This is a complicated expression which becomes very simple at £ = 1/6. In fact 

1^ = ^ = 0. (88) 



This important equation has root in the conformal invariance of the action (I15|) 
at £ = 1/6, as pointed out after eq. (fT5|) . 
Because of ea. ([55|) . we have 

- p<t> + 3p0 = =*> p = ip . (89) 

So one can add the classical contribution of action (|15[) to Friedmann equations 
by considering it as a barotropic fluid with a;^ = 1/3! But we know that adding 
any barotropic fluid with to > — 1 does not change the structure of phase-space, 
so this is also true for Lagrangian (Ti"5j) . This completes our proof. 
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5 Conclusion 



In this paper we study the contribution of quantum phenomena on the structure 
of R — H phase-space of F(R, G) gravity models. It is shown that, both the 
location of de-sitter attractors and their stability properties, can change because 
of quantum corrections. These corrections are not, in general, from the type of 
small, and therefore ignorable, corrections, instead the topology of attractors 
may change. For the case where there is no attractor, it may produce some 
attractors, an example of it is discussed in section 3.1. In some cases, it in- 
creases the number of the attractors and their stability characteristics, with an 
important property: The increased attractors do not converge to classical ones 
at h — ► limit, and more, their locations blow up. A fact which is known in 
the context of " singular perturbation" . An example of this kind of behaviors is 
given in section 3.2. Also it is possible that the quantum contributions change 
the attractors, such that not only they do not converge to classical locations, 
but also the quantum and classical attractors go more distant from each other 
at h — > limit. 

It is also shown that the barotropic perfect fluids with u> > —1, do not 
change our results. Interestingly, the coupled-quintessence Lagrangian, which 
is responsible for producing the quantum terms, is from this type. 

Besides the attractors studied in this paper, there are two other classes of 
attractors in F(R,G) gravity models [14]. One of them is the singularities 
of F(R, G) Lagrangian, which is always lead to the stable attractors, and the 
other class is one known as critical curve. The curve R — 12H 2 is the location of 
infinite number of the stable attractors for the models whose Lagrangians satisfy 
eq. ([3TJj) (with b' = 0). That is, eq.flHJ) (with b' = 0) satisfies for all R and Hs. 
F(R, G) = R 2 g(G/R 2 ),F{R, G) = aG + f3R 2 , F(R, G) = R - 6G/R, • • • , are 
examples of these Lagrangians [14]. Now it can be easily shown that the, so 
called, singular attractors do not change because of the quantum terms pa and 
Pa, but the critical curve R = 12H 2 , reduced to a single point. The reason is 
easy. If F(R, G) satisfies 

^RF R + GF G -F = (90) 

for all R, then for these actions, the equation of critical point (|30)) becomes 

= = ^b'R 2 , (91) 

which has a unique solution R c = 0. So the quantum effects reduce the infinite 
number of stable attractor to an attractor at R c = 0. 

Of course, the critical curve R = 12H 2 still exists in our case, but for 
F(R, G)s which satisfy 

±RF R + GF g -F= h>'R 2 (92) 

for all R and Hs. For these functions, since ea.(|30p satisfies by F(R,G), this 
equation does not impose any constraint on R and H, and ea. (|3ip is the only 
equation which specifies the critical points. So R — 12H 2 becomes the critical 
curve of these models. An example of these kinds of models is 

f(R) = (C+^-\nR)R 2 , (93) 
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which satisfies eg . (j9"2")l . At b' — ► limit, it gives f(R) — R 2 , which is stated 
before (the F(R, G) = aG + (3R 2 model, with a = 0). 
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